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The aim of these notes is both to review the standard understanding of the 
^ ' Hawking effect, and to discuss the modifications to this understanding that might 

\ be required by new physics at short distances. The fundamentals of the Unruh 

effect are reviewed, and then the Hawking effect is explained as a "gravitational 
Unruh effect" , with particular attention to the state-dependence of this picture. 
The order of magnitude of deviations from the thermal spectrum of Hawking 
l/^ \ radiation is estimated under various hypotheses on physics at short distances. 

On ' The behavior of black hole radiation in a linear model with altered short distance 

physics — the Unruh model — is discussed in detail. 

1. Introduction 



X 



When a black hole forms in a collapse process, virtually all features of the original 
collapsing object are erased. All that survives of its origin is the total mass, angu- 
lar momentum, and charges. What remains is a pure object, made of nothing but 
empty spacetime. Thus a non-rotating uncharged black hole is characterized by its 
Schwarzschild radius Rs = 2GM/c^. The length Rs provides the only scale. 

If we bring in h something new happens. The "empty space" that the black hole is 
made of wakes up and is full of vacuum fluctuations of all fields (including, presumably, 
the gravitational field itself). And these vacuum fluctuations are propagating in the 
black hole background. Hawking discovered^ in 1974 that this black hole vacuum is 
not stable. He showed that a thermal flux of positive energy at temperature Th = 
hc/AnRs flows away from the vicinity of the black hole. (For a stellar mass black 
hole this temperature is very low, since the Schwarzschild radius is about a kilometer 
{Th — 10~''K (1.5km/ Rs)), but the temperature could be very high for a much smaller 
black hole that might have formed in the early universe.) Moreover, negative energy 
flows across the horizon into the hole, so total energy in the fixed, static, background 
is conserved. Presumably the gravitational back-reaction to this process diminishes 
the mass of the hole, and the process continues until the hole has lost all of its mass. 



Fig. 1. An outgoing mode typical of the Hawking radiation. Th corresponds to a wavelength of 



The fact that the Hawking radiation has a universal, thermal, spectrum can be 
understood from the observation that there is only one scale in the problem, the 
Schwarzschild radius Rs- Other than indicating the absorption cross-section of the 
black hole, the radiation can carry no information. That is, it has maximum entropy, 
and so must have a thermal spectrum. The temperature corresponds to a wavelength 
of Stt'^Rs- Since h is involved, one might think there is a second scale in the problem, 
the Planck length Lp = {hG/c^f/'^ ~ lO^^^ 

cm. However, G does not enter into 
Hawking's result because the gravitational interaction plays no role. Only when one 
includes the back-reaction does G or Lp enter. Indeed, the lifetime of the hole is of 
order R\/ L'\,c. (This is roughly 10^^ times the present age of the universe for a stellar 
mass black hole. For a primordial black hole of mass ~ lO^^gm however the lifetime 
is only the present age of the universe.) 

Upon further reflection however, it is not so clear that the Planck scale physics 
is irrelevant to black hole radiationiffli. Consider a vacuum fluctuation mode that 
ends up at inflnity populated with Hawking radiation. This mode emerges from the 
vicinity of the horizon, and must climb out of a very deep well (see Fig. |l|). It 
is tremendously redshifted on the way out. If the mode has a wavelength \out when 
observed a time t after the hole formed, then its initial wavelength Aj„ must be of order 
exp(— t/2i?s)Ao„(, which is very quickly way below the Planck length. For instance 
if Xout ~ Rsi then Aj„ ~ Lp when t ~ 2Rsln{Rs/Lp), which is about a millisecond 
for a solar mass black hole. Hawking's original calculation assumed that vacuum 
fluctuations are scale invariant down to such arbitrarily sub-Planckian dimensions. If 
this assumption is no good, then the true behavior of quantum flelds around a black 
hole could be quite different. 

The aim of these notes is both to review the standard understanding of the Hawk- 
ing effect, and to discuss the modiflcations to this understanding that might be re- 
quired by new physics at short distances. In particular, I flrst review the Unruh 
effect, namely, the fact that an accelerated observer in Minkowski spacetime sees the 
vacuum as a thermal state. This is then used to give an elementary explanation of 
the origin of the Hawking effect, and to argue that, whatever short distance fleld dy- 



namics occurs on sub-Planckian scales, it will give rise to the Hawking effect as long 
as it produces an ordinary vacuum structure for the outgoing modes near the horizon 
at scales A satisfying Lp < \ <ti Rs- The order of magnitude of deviations from the 
thermal spectrum of Hawking radiation is estimated under various hypotheses. In the 
remainder, the behavior of a linear model with altered short distance physics — the 
Unruh model — is discussed. This model consists of a scalar field propagating on a 
black hole background with non-Lorentz-invariant higher derivatives in the action. 
The higher derivatives lead to a nonlinear dispersion relation at large wavevectors, 
but the field equations remain linear. 

In the following, unless otherwise indicated, units are adopted for which h = c = 
G = kB = l. 



2. The Unruh effect 



A uniformly accelerated observer in fiat spacetime perceives the vacuum state of 
quantum fields as a thermal state with temperature Tu = a/27r, where a is the proper 
acceleration. This is called the Unruh effect!, and Tu is the Unruh temperature. To 
get a feeling for the magnitude of this temperature, we can put h and c back in: 

SO Tf/ is roughly h divided by the time it takes for the velocity to change by the speed 
of light. 

The Unruh effect can be understood as a consequence of the invariance of the 
vacuum under the symmetries of Minkowski spacetime, translations and boosts, to- 
gether with the spectral condition that the energy of physical states be nonnegative. 
I will give two derivations of this effect below, both of which are valid for arbitrary 
interacting scalar fields in spacetime of any dimension. (The generalization to fields of 
nonzero spin is essentially trivial.) The second derivation, which uses a path integral 
method, is much slicker and generalizes to the Hartle-Hawking state for black holes 
and all other curved but static spacetimes with bifurcate Killing horizons. My reasons 
for giving the first derivation as well are threefold. First, it is very concrete, second, 
it illustrates very nicely the dangers of being sloppy with analyticity properties, and 
third, I have not seen it done in the literature in this much generality in this simple 
manner. A brief guide to the literature on these matters is provided in subsection 
below. 

The Minkowski line element in two dimensions can be written in both "Cartesian" 
(Minkowski) and "polar" (Rindler) coordinates: 

ds^ = dt" - dz^ = i'^drf - rfe' (2) 

where the coordinates are related by 

t = ^sinh.rj, z = ^cosh.rj. (3) 



The coordinates (?7, nonsingular in the ranges C, G (0, oo) and r] G (—00,00), 
and cover the "Rindler wedge" x > |t| in Minkowski space (see Fig. 0). In the first 
form of the hne element the translation symmetries generated by the Killing vectors 
d/dt and d/dz are manifest, and in the second form the boost symmetry generated 
by the Killing vector d/drj is manifest. The latter is clearly analogous to rotational 
symmetry in Euclidean space. The full collection of translation and boost symmetries 
of Minkowski spacetime is called the Poincare group. 




Fig. 2. Two-dimensional flat spacetime in Minkowski and Rindler coordinates. A hyperbola of 
constant f is a uniformly accelerated timelike worldline with proper acceleration A boost shifts 
77 and preserves ^. 



2.1. Two-point function and KMS condition 

A thermal density matrix p = exp{—j3H) has two identifying properties: 
First, it is obviously stationary, since it commutes with the Hamiltonian H. Second, 
because exp(— coincides with the evolution operator exp{—itH) for t = —iP, 
expectation values in the state p possess a certain symmetry under translation by 
—if] called the KMS conditioni'i: Let (A)^ denote the expectation value tr{pA), and 
let At denote the time translation by t of the operator A. Using cyclicity of the trace 
we have 

{A^^,B), = Z-Hr{e-^^ie^^Ae-^'')B) (4) 
= Z-Hr{e-^"BA) (5) 
= {BA)p. (6) 

Note that for nice enough operators A and B, {A_iT-B)p will be analytic in the strip 
< r < /3. Now let us compare this behavior with that of the two-point function 
along a uniformly accelerated worldline in the Minkowski vacuum. 



If, as is usual, the vacuum state shares the symmetry of Minkowski spacetime, 
then, in particular, the 2-point function G{x,x') = {(j){x)(f){x')) must be a Poincare 
invariant function of x and x' . Thus it must depend on them only through the 
invariant interval, so one has G{x^ x') = f{{x—x'Y) for some function /. Now consider 
an "observer" traveling along the hyperbolic trajectory ^ = a~^. This worldline has 
constant proper acceleration a, and at] is the proper time along the world line. Let 
us examine the 2-point function G{ri,r]') = G{x{ri),x{ri')) along this hyperbola. It is 
clearly stationary, due to the boost invariance of the vacuum, so it can only depend 
on 1] and rj' through the difference rj — rj'. Therefore 

G(r?,r?') = G{x{v-v'),xm (7) 

= fiixiv-v')'mr) (8) 
= /(4a-2sinh2[(r/-r^')/2]), (9) 

where the third equality follows from @. Now, since sinh^(?7/2) is periodic under 
translations of t] by 27Ti, it appears that G{r], rj') is periodic under such translations in 
each argument. In terms of the 2-point function the KMS condition implies G{ri,ri' + 
i(3) = G{ri',ri), which is not the same as translation invariance in each argument. 
Does this mean that in fact the 2-point function in the Minkowski vacuum along the 
accelerated worldline is not thermal? The answer is "no" , because the above "proof 
that G(77, t]') is periodic was bogus. First of all, a Poincare invariant function of x and 
x' need not depend only on the invariant interval. It can also depend on the invariant 
step-function 6'(x° — x"^)6{{x — x')"^). More generally, the analytic properties of the 
function / have not been specified, so one cannot conclude from the periodicity of 
sinh^(r7/2) that / itself is periodic. For example, / might involve the square root, 
sinh(?7/2), which is anti-periodic. In fact, this is just what happens. 

To reveal the analytic behavior of G{x,x'), it is necessary to incorporate the 
condition that the spacetime momenta of states in the Hilbert space lie inside or on 
the future light cone. One can show (by inserting a complete set of states between 
the operators) that this implies there exists an integral representation for the 2-point 
function of the form 

G{x,x') = J d'''ke{k°)J{k^)e-'''^'-'-^'\ (10) 

where J{k'^) is a function of the invariant k"^ that vanishes when k is spacelike. Now 
let us evaluate G{ri, rj') along the hyperbolic trajectory. Thanks to boost invariance, 
it suffices to consider rj' = 0. Thus we have 

G(r/,0) = J rf"A;0(A;°)J(A;2)e-*''''['="'^^°'^'^-^'(™^^''-^)l. (11) 

Now it still looks as though this is periodic in translation of rj by i2TTm for any integer 
m, but closer inspection reveals the difficulty that the integral fails to be convergent 



if the imaginary part of sinh rj becomes positive. One has 



sinh(?7 — 2^^) = sinh r/ cos 6* — 2 cosh 7] sin ^ (12) 
cosh (r] — 20) = cosh cos — i sinh ?7 sin 0. (13) 

Since > \k^\ the integral is seen to converge for < < vr, but one cannot go 
past IT. To see if the function has an analytic extension beyond 6 = it and, if so, 
to determine the value at 6 = 27r, one needs a different representation of G{ri). A 
suitable representation can be obtained by observing that Lorentz invariance allows 
us to transform to the frame in which x — x' has only a time component which is 
given by the invariant norm [{x — x'^Y^'^ = 2a^^ sinh(r7/2). Thus we have 

G{t], 0) = y 0(fcO) J(fc2)g-i2a-lfc0sinh(r,/2)_ ^^4) 

In this form it is clear that r] — iO can be extended all the way down to 9 = 2tt. 
Since sinh[(r7 — z27r)/2] = — sinh(77/2) = sinh(— 77/2), we can finally conclude that 
G{r] - i2n, 0) = G{-r], 0) = G(0, rj), which is the KMS condition (§). 

2.2. The vacuum state as a thermal density matrix 

The essence of the Unruh effect is the fact that the density matrix describing the 
Minkowski vacuum, traced over the states in the region 2; < 0, is precisely a Gibbs 
state for the boost Hamiltonian Hb at a "temperature" T = l/27r: 

Tr,<o|0)(0| = Z-'exp{-2nHB), (15) 

Hb = I Tatid/drirdE' (16) 

This rather amazing fact has been proved in varying degrees of rigor by many different 
authors. A sloppy path integral argument making it very plausible will be sketched 
below. 

Since the boost Hamiltonian has dimensions of action rather than energy, so does 
the "temperature" . To determine the local temperature seen by an observer following 
a given orbit of the Killing field, note from (0) that the norm of the Killing field d/dr] 
on the orbit ^ = is a~^, whereas the observer has unit 4- velocity. If the Killing 
field is scaled by a so as to agree with the unit 4- velocity at ^ = a~^, then the boost 
Hamiltonian (p!6D and temperature are scaled in the same way. Thus the temperature 
appropriate to the observer at ^ = a^^ is T = a/27r. Since a is the proper acceleration 
of this observer, we recover the Unruh temperature defined above. Alternatively, the 
two-point function defined by (|15|) along the hyperbola obviously satisfies the KMS 
condition relative to boost time t] at temperature l/27r. When expressed in terms of 
proper time arj, this corresponds to the temperature a/27r. 



One can view the relative coolness of the state at larger values of ^ as being due 
to a redshift effect — in this case a Doppler shift — as follows. Suppose a uniformly 
accelerated observer at sends some of the thermal radiation he sees to another 
uniformly accelerated observer at > ^o- This radiation will suffer a redshift given 
by the ratio of the norms of the Killing field: say p is the spacetime momentum of 
the radiation. Then p ■ {d/dif) is conserved!, but the energy locally measured by the 
uniformly accelerated observer isp- {d / drj) /\d / drj\, so that Ei/Eq = \d/dri\o/\d/dri\i. 
This is precisely the same as the ratio Ti/Tq of the locally measured temperatures. 
At infinity \d/dri\ = C, diverges, so the temperature drops to zero, which is consistent 
with the vanishing acceleration of the boost orbits at infinity. 

The path integral argument to establish (|l^) goes like this: Let H be the Hamil- 
tonian generating ordinary time translation in Minkowski space. The vacuum |0) is 
the lowest energy state, and we suppose it has vanishing energy: H\0) = 0. If \ip) is 
any state with nonzero overlap with the vacuum, then exp{—TH)\ip) becomes propor- 
tional to |0) as r goes to infinity. That is, the vacuum wavefunctional \l/o[0] for a field 
(f) is proportional to exp{—TH)\il)) as r — *• oo. Now this is just a matrix element of 
the evolution operator between imaginary times r = — oo and r = 0, and such matrix 
elements can be expressed as a path integral in the "lower half of Euclidean space: 



^o[0]= / P0exp(-J) (17) 
where / is the Euclidean action. 



The key idea in recovering ([15D is to look at ([T7|) in terms of the angular "time"- 



slicing of Euclidean space instead of the constant r slicing. (See Fig. |^.) The relevant 



Fig. 3. Time slicings of Euclideanized Minkowski space. The horizontal lines are constant r surfaces 
and the radial lines are constant 9 surfaces. 

Euclidean metric (restricted to two dimensions for notational convenience) is given 
by 

ds^ = dr^ + da^ = p^dO"^ + dp^. (18) 



Adopting the angular slicing, the path integral ( p!7D is seen to yield an expression for 
the vacuum wavefunctional as a matrix element of the boost Hamiltonian (|T6|) which 
coincides with the generator of rotations in Euclidean space: 



(0L0ii|O) = Ar(0i| exp(-7r/7B)|0ij), (19) 

where (pL and 0r are the restrictions of the boundary value 0(0) to the left and right 
half spaces respectively, and a normalization factor M is included. The Hilbert space 
Ti^ on which the boost Hamiltonian acts consists of the field configurations on the 
right half space ^ > 0, and is being identified via reflection (really, by reflection 
composed with CPtEI'EI) with the Hilbert space TCl of field configurations on the left 
half space z < 0. The entire Hilbert space is H = Hl ® Ti-R, modulo the degrees of 
freedom at 2; = 0. (The boundary conditions at 2; = are being completely glossed 
over here.) Now consider the vacuum expectation value of an operator that is 
localized on the right half space: 

(0|(9h|0) = (O|0^0«)(0^0«|O«|(/.'^0^)(</.i^</.'«|O) (20) 
= Af^M expi-7THB)\(f)R){(f)R\OR\(f)'j,){(j)'j,\ expi-TTHsML) (21) 
= AfHr{e"^^''^OR), (22) 



where summation over intermediate states is implicit, and (19) was used in the second 
equality. This shows that, as far as observables located on the right half space are 
concerned, the vacuum state is given by the thermal density matrix (P!3|). More 
generally, this holds for observables localized anywhere in the Rindler wedge. 

This path integral argument directly generalizes to all static spacetimes with a 
bifurcate Killing horizon, such as the Schwarzschild and deSitter spacetimeJii'0. In 
the general setting, the state defined by the path integral cannot be called "the" 
vacuum, but it is a natural state that is invariant under the static Killing symmetry 
of the background and is nonsingular on the time slice where the boundary values of 
the field are specified, including bifurcation surface. 

2.3. Some references 

The fact that the Minkowski vacuum is a thermal state for the boost Hamilto- 
nian was proved in axiomatic quantum field theory by Bisognano and WichmannEl, 
as a theorem about the action of complex Lorentz transformations on the vacuum. 
The relevance of this theorem to the Unruh and Hawking effects was recognized by 
Sewelli, who generalized the framework to curved spacetimes. In completely inde- 
pendent work (as far as I know) the path integral argument has been given by many 
authors, perhaps the first beingJJnruh and Weis^^. The review articles by TakagicJ 
and by Fulling and Ruijsenaarsc3 cover various aspects of the relation between accel- 



eration and temperature in quantum field theory, and contain many other references. 

3. Hawking radiation 

The Hawking effect can be understood as a kind of "gravitational Unruh effect" 
as followsi'0. Consider a static, spherically symmetric, asymptotically flat black hole 
spacetime. This spacetime has a time translation Killing field x which is a unit 
timelike vector at infinity, becomes null on the horizon, and is spacelike within the 
horizon. An observer at fixed radius is uniformly accelerated, with an acceleration 
that vanishes at infinity and diverges as the horizon is approached. What can be said 
about the state of the outgoing modes as viewed from infinity in such a spacetime? It 
should be consistent with redshifting from the state of the same modes viewed from 
a vantage point closer in to the black hole, however that of course does not determine 
the state. 

In order to infer the existence of Hawking radiation one must assume some bound- 
ary condition on the state of the quantum field. In principle, this should just be the 
assumption that the initial state was not too pathological and contained matter that 
would later collapse and form a black hole. However, as discussed in the introduction, 
to connect this initial condition to the final state one must evolve initially arbitrarily 
short wavelength field modes through the collapse. Instead let us first consider a 
boundary condition that does not involve us in the physics of the trans-Planckian 
domain. Thus let us suppose that the state just outside the horizon, after the black 
hole has formed, looks like the Minkowski vacuum at very short, but longer than 
Planckian, distances. If this is the case, then a static observer hovering just outside 
the horizon with a tremendous acceleration should experience the Unruh effect, and 
will perceive the state to be thermal at some very high temperature. This thermal 
state is then redshifted to infinity, as previously discussed in the case of fiat space- 
time. However a crucial difference now arises: the norm of the Killing field x is 
finite at infinity, rather than divergent. Thus the Unruh radiation survives its trip 
to infinity, where it arrives with a non-zero temperature. Let us see how this works 
quantitatively. 

For concreteness let us work with a Schwarzschild black hole, whose line element 
is given in Schwarzschild coordinates by 

ds^ = (1 - 2M/r)dt^ - (1 - 2M/r)-^dr^ - r\de^ + sin^ed(j)'^). (23) 
The static Killing field is x = d/dt, and its norm is given by 

iV(r) = (x-x)'/' = (l-2M/r)V2. 

is also called the "lapse" function, since it determines how much proper time 
ds corresponds to a given coordinate time dt via ds = Ndt at fixed {r,9,(j)). The 



acceleration of a worldline at constant radius and angles is found to be 



a = N-^{1- N^K, 



(24) 



where k = 1/4M is the "surface gravity" of the hole. At the horizon the acceleration 
diverges and at infinity it vanishes. The Unruh temperature corresponding to the 
acceleration is a/2'ir. As viewed from infinity this temperature suffers a redshift given 
(as explained in subsection ) by the ratio N{r)/N{oo) = N{r). Thus the Unruh 
temperature for an observer at r, as viewed from infinity, is 



where Th = k/27t = l/SvrM is the Hawking temperature. As r approaches the horizon 
the redshifted Unruh temperature approaches the Hawking temperature, while as r 
goes to infinity it vanishes. 

Since the value of Tu^oo{r) in (|25|) very much depends on r, we appear to have 
an inconsistency. How can the state viewed from infinity have many different tem- 
peratures? It cannot, of course. The inconsistency arises only if we assume that at 
each radius the static observer experiences the usual Unruh effect, but this assump- 
tion requires that along each of these worldlines the 2-point function of the quantum 
field has the same form as it would have in Minkowski space. We have no basis for 
making such an assumption, except in the limit of extremely short distances where 
we expect all states to coincide. Since a higher Unruh temperature arises from a 
shorter distance probe, this assumption should get better and better as the static 
world line approaches the horizon. In this limit, Tu^oo{r) approaches the Hawking 
temperature, so it is the Hawking temperature that is selected by the Minkowskian 
boundary condition on the quantum state at very short distances. 

There are of course many other states of the quantum field. For example, there 
is the "Boulware vacuum", i.e., the Fock vacuum devoid of positive Killing-frequency 
excitations. In this state the static observers see no particles at any radius. However, 
this state is not generally believed to be the state that arises when a black hole forms 
in a collapse process. In fact if one evolves the state of a free quantum field on the 
background of a collapsing black hole, one sees that the final state is the "Unruh 
vacuum". Near the horizon these two states are tremendously different. Whereas 
the Unruh vacuum looks like the Minkowski vacuum at very short distances, the 
Boulware vacuum does not. For example, the renormalized expectation value of the 
stress-energy tensor in the Boulware vacuum blows up as the horizon is approached. 

4. Imposing a cutoff 

Now imagine there is a short distance cutoff of some kind on the validity of 
ordinary field theory, so one cannot impose the vacuum condition as above in the 
limit of arbitrarily short distances. What effect would that have on the Hawking 



(25) 



effect? Of course the answer depends on what kind of short distance physics is 
imagined. Let us suppose first that, whatever the details at the cutoff, the vacuum 
boundary condition is still satisfied to a fair degree at scales longer than the cutoff Ic 
but much shorter than the Schwarzschild radius. That is, we assume the out vacuum 
at scales A satisfying 

< A < Rs. (26) 

One would then expect the Hawking effect to still occur, but with small deviations 
from the usual Hawking spectrum. The presence of the cutoff presumably limits 
how close to the horizon the accelerated observer can be assumed to experience the 
Unruh effect. If the minimum allowed lapse is A^min, then we might expect from 
(|25| ) deviations at least of order iV^j^ from the usual Hawking fluJ, since even if the 
observer at A^min experiences the usual Unruh effect (which is presumably not quite 
the case), the redshifted temperature ( |2^ ) will not agree with Th- The value of A^min 
would depend on the nature of the cutoff. 

4.I. Invariant cutoff 

Without breaking local Lorentz invariance, one can imagine a cutoff that limits 
the acceleration for which the usual Unruh effect takes place. This might be the case 
if there is a minimum proper time that should be considered along the worldline. If 
the acceleration (^) is required to be less than then A^min ~ ^^^c — h/Rs- In 
this case the above argument suggests deviations of order {1^/ Rg)^ from the usual 
Hawking effect. This might be relevant, for instance, in string theory. 

4-2. non-Lorentz-invariant cutoff 

Lorentz invariance of the laws of physics has been confirmed to some precision up 
to some finite boost factor 7 relative to the cosmic rest frame. The boosts that play 
a role in the Hawking effect however grow exponentially without bound in time as 
7 ~ exp{t/2Rs). Normally physicists do not have confidence extrapolating infinitely 
beyond their observations, but the case of Lorentz invariance seems to be an exception. 
Allowing this exception seems to force us into assuming that there are an infinite 
number of field modes in any finite spatial region, no matter how small. And all of 
these modes play a role in the Hawking effect as usually formulated. It is therefore 
worthwhile to ask what would happen to the Hawking radiation in the presence of a 
non-Lorentz-invariant cutoffBi. This question will be considered in a detailed model 
later in these notes. For now, let us just assume as above that the physics produces 
the out vacuum at scales A satisfying (^) in the free-fall frame of the black hole, that 
is, in the frame carried in to the black hole by observers who fall freely from far away, 
starting at rest. 

Of course this definition of the preferred frame is ambiguous. In the spherically 



Fig. 4. The null cutoff wave-vectors and the four-velocities of the freely falling and static observers. 
The diagram is drawn in the free-fall frame, so the static observer appears to be moving outward at 
high speed. 



symmetric case we can restrict to radial trajectories, and consider only those trajecto- 
ries that are always infalling, ignoring those that pass through the collapsing matter 
and later fall back. More generally, the fundamental cutoff theory would be needed 
to even formulate the nature of the cutoff. But in order to explore the qualitative 
implications of such a cutoff, it seems reasonable to stick to the spherically symmetric 
case and to impose the cutoff in the falling frame defined above. 

The cutoff for the static observer is Doppler shifted relative to the cutoff for the 
free-fall observer. The outgoing modes are redshifted while the ingoing modes are 
blueshifted. Let ^cin and /cc,out denote ingoing and outgoing null wave-vectors with 
the cutoff wavelength, and let us and -Ugtat denote the unit 4-velocities of the free- 
fall and static observers (see Fig. ^). Then the cutoff frequency is given in the 
free-fall frame by ujc = ^c,in ■ = ^cout ■ u^, while it is given in the static frame by 
^c!kf = ^c,in ■ Mstat and cu^^^ut = ^c,out " ^stat for the in- and out-going modes respectively. 
Using these definitions near the horizon where -C 1 one finds 



For the Unruh effect we must require that the acceleration timescale a is longer 
than the cutoff in the accelerated frame. For the outgoing modes, which are relevant 
for the Hawking effect, this means that N/k > N~^lc. This non-Lorentz-invariant 
cutoff thus leads to a much larger minimum lapse, A^min — ih/RsY^'^- In this case 
one might expect deviations from the Hawking flux of order Ic/Rs, which is still very 
small for all but Planckian holes. 

A quantum field theoretic calculation of the out state at infinity starting from 
the out vacuum near the horizon at scales (p6D in the free-fall frame has also been 
carried outi. To do this calculation it is necessary to work with localized wavepackets. 
Being careful about the errors caused by the inevitable "tails" of the wavepackets, 
an upper bound of order {lp/Rs)~^^'^ for the deviations from the Hawking spectrum 
was found, which is much larger than the estimate obtained from the Unruh effect 




(27) 
(28) 



argument above. However, this was only an upper bound. 

4.3. Horizon fluctuations 

So far I have been speaking of a cutoff in terms of a modification of the field theory 
on a fixed background. On the other hand, taking into account the back reaction and 
quantum fluctuations in geometry one expects the horizon to fluctuate in some sense. 
It is interesting to try to estimate the size of these fluctuations and their possible effect 
on the Hawking radiation. Let us consider the fluctuations in the horizon that would 
be expected if the horizon is viewed as a system of Planck areas fluctuating about an 
equilibrium configuration. In equilibrium the typical entropy fiuctuation is given by 
5S ^ 1. Since S — A/A (where A is the horizon area), we therefore expect 5^4 ~ 1. 
Now if we assume the horizon consists of A independent Planck areas with random 
area fluctuation 6a ~ R6r/A 6r/R (in four-dimensional spacctime), then by the 
law of large numbers we have 5A ~ A^^'^5a ~ 5r, so (5r ~ 1 = Ip. This reasoning 
suggests that the horizon is fuzzy on scales of order Ip in the radial coordinate. If 
this is interpreted to mean that one should not apply the acceleration temperature 
argument any closer to the horizon than r = Rs + lp, one finds for the minimum lapse 
A^min ~ {Ip/RsY''^. This is the same result as one would have obtained by imposing 
a cutoff at /p in the free-fall frame as above. 

5. Models with a non-Lorentz-invariant cutoff 

In the previous section we argued crudely that the Hawking effect will occur 
independently of short distance physics as long as the out vacuum boundary condition 
holds near the horizon for wavelengths satisfying /c < A <^ -Rs, up to small corrections 
for large black holes. This leaves us with the following questions. Given a particular 
theory with some new physics at short distances, 

1. Does the abovementioned out vacuum boundary condition hold? 

2. Exactly how large are the deviations from the the thermal Hawking fiux? 

3. Are the deviations from the thermal Hawking flux small even at very short 
wavelengths? 

4. Do the deviations for short wavelengths accumulate to make a large difference 
in any physical quantity, such as the energy flux or energy density? 

A class of models in which these questions can be quantitatively addressed with 
relatively simple methods is obtained if the usual wave equation for the quantum 
fleld is modified by the addition of higher derivative terms, while linearity of the 
fleld equation is preserved. Such a modification can be covariant, or it can be done 



in a manner that violates local lorentz invariance. I know of no suitable covariant 
model for ordinary fields, since covariant higher derivatives always lead to negative 
norm ghosts. However, string theories can be Lorentz-covariant and ghost-free, and 
it might be feasible and interesting to study Hawking radiation in non-interacting 
string theory. Attempts in this direction have been made but so far none have been 
able to compute the spectrum of black hole radiation. On the other hand, several 
non-covariant models ones have been studied recently. The mother of them all was 
Unruh's sonic black hole. Although the subsequent models can all be interpreted in 
terms of ordinary black holes without reference to fluid flow, I will flrst describe the 
sonic black hole model, since it is very interesting and provides good motivation for 
the other models. 



5.1. Sonic black hole 

Unruhi pointed out in 1981 that perturbations of an irrotational, inviscid, barotropic 
(p = fluid flow behave like a massless fleld propagating in a curved spacetime 

whose metric is determined by the background flow. If the flow is supersonic some- 
where, then there is a horizon in the effective spacetime geometry, and one has a model 
of a black hole. Quantizing the fluid perturbations, Unruh argued that, provided the 
out vacuum boundary condition holds just outside the sonic horizon, the sonic black 
hole will radiate thermal phonons at the temperature Ks /27r, where Kg is the gradient 
of the velocity fleld at the horizon. (Numerically, Ks/2ti ~ 10~^K (|^/^^~~)-) This 
is already interesting, but the real reason Unruh invented this model was in the hope 
of studying both the consequences of the molecular nature of the fluid and the quan- 
tum backreaction. Since the molecular structure has a preferred local frame — the 
comoving frame of the fluid — the molecular fluid model lacks covariance in the fluid 
metric. That is, the molecular fluid model displays a non-Lorentz-invariant cutoff. 

Consider the sonic black hole with liquid Helium-4 as the fluidS. The sound 
fleld describes the quasiparticle excitations, and satisfles the usual wave equation 
for long wavelengths. The atomic nature of the fluid begins to show itself in the 
nonlinearity of the quasiparticle dispersion relation uj{k) at shorter wavelengths. This 
nonlinearity implies that the group velocity is fc-dependent, at flrst dropping as k 
grows. This has important implications for the propagation of wavepackets in the 
sonic black hole background. In particular, one can infer that an outgoing wavepacket 
propagated backwards in time toward the sonic horizon will be blueshifted, hence its 
comoving group velocity will drop. The packet will not quite make it all the way to 
the horizon, but will stop where its comoving group velocity is equal to the negative 
of the background flow velocity. This is an unstable situation, and what happens at 
this stage is remarkableil: the packet is blueshifted a bit more, the comoving group 
velocity drops a bit more, and the packet is swept back away from the horizon by 
the (time-reversed) flow! This is very different from the behavior of a wavepacket 



satisfying the ordinary wave equation, which is to get squeezed more and more up 
against the black hole horizon, while exponentially blueshifting. 

5.2. Unruh model 

It is not necessary to deal with liquid Helium to study the behavior just described. 
All that is essential is the curvature of the dispersion relation, and this can be modeled 
in linear field theory with higher spatial derivative terms0. Since in the case of a fluid 
the dispersion relation is physically specified in the comoving frame, the comoving 
time derivative should be left unchanged. Unruh studied the Hawking process in a 
model of this sort, by numerically integrating the partial differential equation govern- 
ing the propagation of a wavepacket . Because of the reversal of wavepacket motion at 
the horizon, it is possible to deduce the final state from an initial condition on ingoing 
modes far from the horizon, even in a stationary background (see Section for further 
discussion of this issue). Unruh found, to numerical accuracy of his computation, that 
the usual Hawking flux is emitted from the hole. Subsequently the same and similar 
models were studied using using both analytical approximationJ^l, and by first using 
stationarity to reduce the problem to the solution of ordinary differential equationfll. 
Using the last described method, it was possible to reliably compute deviations from 
the usual thermal spectrum of Hawking radiation. 

Here I will describe the Unruh model and some of its generalizations without 
further reference to the fluid flow analogy. The model consists of a free, hermitian 
scalar field propagating in a two dimensional black hole spacetime. The dispersion 
relation for the field lacks Lorentz invariance, and is specified in the free fall frame 
of the black hole, that is, the frame carried in from the rest frame at infinity by 
freely falling trajectories. Let denote the unit vector field tangent to the infalling 
worldlines, and let s" denote the orthogonal, outward pointing, unit vector, so that 
gai3 _ y^a^iB _ ^0^/3 ^g^g pjg_ ^ rjj^g actiou is assumed to have the form: 

S = ^Jd'x ^gg''^V^ct>*Vpct>, (29) 
where the modified differential operator is defined by 

= m"(9„ (30) 
= F(s"a«). (31) 

The time derivatives in the local free fall frame are thus left unchanged, but the 
orthogonal spatial derivatives are replaced by F{s°'da)- The function F determines 
the dispersion relation. For the moment it will be left unspecified. Invariance of 
the action (|29|) under constant phase transformations of guarantees that there is a 
conserved current for solutions and a conserved "inner product" for pairs of solutions 
to the equations of motion. However, since T>a is not in general a derivation, simple 
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Fig. 5. A patch of spacetime showing a free-fall trajectory and some t and x (Lemaitre-like) coordi- 
nate lines, u and s are orthonormal vectors, and the derivative along s is modified, while that along 
u is just the partial derivative. The notations St and dx denote d/dt and d/dx respectively, and St 
is the Killing vector. 



integration by parts is not allowed in obtaining the equations of motion or the form 
of the current. We shall obtain these below after further specifying the model. 
The black hole line elements we shall consider are static and have the form 

ds"^ = dt^ - {dx - v{x) dtf. (32) 

This is a generalization of the Lemaitre line element for the Schwarzschild spacetime, 
which is given by v{x) = —\j2M/x (together with the usual angular part). We shall 
assume f < 0, dv/dx > 0, and f — fo as x — > oo. dt is a. Killing vector, of squared 
norm 1 — v'^, and the event horizon is located at v = —1. The curves given by 
dx — vdt = are timelike free fall worldlines which are at rest (tangent to the Killing 
vector) where v = 0. Since we assume v < these are ingoing trajectories, v is their 
coordinate velocity, t measures proper time along them, and they are everywhere 
orthogonal to the constant t surfaces. (See Fig. ^.) We shall refer to the function 
v{x) as the free-fall velocity. The asymptotically flat region corresponds to x — > oo. 

In terms of the notation above, the orthonormal basis vectors adapted to the free 
fall frame are given hj u = dt + vdx and s = dx, and and in these coordinates g = —1. 
Thus the action (^Of ) becomes 

S = ^J dtdx (\{dt + vdx)(P\^ - \F{dx)<P\^). (33) 

If we further specify that F{dx) is an odd function of d^, then integration by parts 
yields the field equation 



{dt + dxV){dt + vdx)<P = F\dx)<P. 



(34) 



The conserved inner product in this case is given by 

{(j),ij)=t J dx{<j)*{dt + vd.,)^-^{dt + vd^)<f)*), (35) 

where the integral is over a constant t shce and is independent of t if (p, satisfy the 
field equation (^41). The inner product can of course be evaluated on other slices as 
well, but it does not take the same simple form on other slices. 

The dispersion relation for this model in flat spacetime, or in the local free fall 
frame, is given by 

uj^ = F\k), (36) 

where F{k) = —iF{ik). Unruh's choice for the function F{k) has the property that 
F'^{k) = k"^ for k <^ ko and F'^{k) = k^ for k ^ ko, where ko is a wavevector 
characterizing the scale of the new physics. We usually think of ko as being around 
the Planck mass. Specifically, he considered the functions 

F{k) = fco{tanh[(A;/fco)"]}". (37) 

Of course there are many other modifications one could consider. Perhaps the simplest 
is given by 

F\k) = e- k^/kl (38) 

This dispersion relation has the same behaviour for small k as the one above, but 
behaves quite differently for large k. It is the one which was studied in Ref. 0. These 
two dispersion relations are plotted in Fig. ^ along with the dispersion relations for 
the ordinary wave equation and for elementary excitations of liquid helium-4. 




5.2.1. Quantization 



To quantize the field we assume the field operator (j){x) is self-adjoint and satisfies 
the equation of motion ( ^4]) and the canonical commutation relations. In setting up 



the canonical formalism, it is simplest to use the time function and evolution vector 
for which only first order time derivatives appear in the action. (Otherwise one must 
introduce extra momenta which are constrained, and then pass to the reduced phase 
space.) This just means that we define the momenta by 

7r = 5L/5{dt<P) = {dt + vd^)(t), (39) 

i.e., TT is the time derivative along the free-fall world lines. The equal time canonical 
commutation relations are then [0(x),7r(y)] = i6{x,y), as usual. 

We define an annihilation operator corresponding to an initial data set / on a 
surface S by 

a(/) = (/,0), (40) 

where the inner product is evaluated on S. If the data / is extended to a solution 
of the field equation then we can evaluate the inner product in (^) on whichever 
surface we wish. The hermitian adjoint of a(/) is called the creation operator for / 
and it is given by 

a\f) = -{rA). (41) 

The commutation relations between these operators follow from the canonical com- 
mutation relations satisfied by the field operator. The latter are equivalent to 

[a{f),a\g)] = {f,g), (42) 

provided this holds for all choices of / and g. Now it is clear that only if / has 
positive, unit norm are the appelations "annihilation" and "creation" appropriate 
for these operators. From (^) and the definition of the inner product it follows 
identically that we also have the commutation relations 

[a{f),a{g)] = -{f,g*), Wif),a^{g)] = -{r,g). (43) 

A Hilbert space of "one-particle states" can be defined by choosing a decompo- 
sition of the space S of complex initial data sets (or solutions to the field equation) 
into a direct sum of the form S = SpQ) Sp*, where all the data sets in Sp have positive 
norm and the space Sp is orthogonal to its conjugate Sp*. Then all of the annihilation 
operators for elements of Sp commute with each other, as do the creation operators. 
A "vacuum" state |\E') corresponding to Sp is defined by the condition a(/)|\E') = for 
all f in Sp, and a Fock space of multiparticle states is built up by repeated application 
of the creation operators to |\E'). 

It is not necessary to construct a specific Fock space in order to study the physics of 
this system. In fact, any individual positive norm solution p defines annihilation and 
creation operators and a number operator N{p) = a^(p)a(p). The physical significance 
of the number operator depends of course on the nature of p. 

There are two types of positive norm wavepackets in which we are interested. 
The first are those corresponding to the quanta of Hawking radiation. These have 



positive Killing frequency, that is, they are sums of solutions satisfying dtcj) = —iuj(j) 
with uj > 0. It is not obvious that such solutions have positive norm in the inner 
product (p5D, and in fact they do not in general. However, using the fact that the 
Killing frequency is conserved, we know that if a positive Killing frequency wavepacket 
were to propagate out to infinity (or any other region where v = 0), the integral for 
its norm would be manifestly positive. Since the norm is conserved, this suffices. 

The other type of positive norm wavepackets we shall employ are those which 
correspond to particles as defined by the free-fall observers. These have positive free- 
fall frequency, that is, they are sums of solutions satisfying {dt + vdx)(p = —ito'cj), 
with u' > 0, on some time slice. These have manifestly positive norm, although the 
free-fall frequency is not conserved. 



5.2.2. Computing the particle creation rate 

In this subsection we obtain the explicit expression for the particle creation rate 
in terms of the norm of the negative frequency part of the ingoing wavepacket that 
corresponds to a given outgoing wavepacket. Let ipout denote a wavepacket solution 



of the field equation (|3^) which at late times is outgoing and localized in the constant 
velocity region, where it has only positive free-fall frequency components, and only 
positive Killing frequency components. Propagating this data backwards in time, 
a part will refiect off the geometry outside the black hole and arrive back in the 
constant velocity region again as a wavepacket ipr with wavevector components at 
the small negative k{uj) root with positive free-fall frequency. (If the ordinary two- 
dimensional wave equation were satisfied, this reflected part would vanish due to 
conformal invariance. However the field equation (|3^) lacks conformal invariance at 
high wavevectors, so in general there will be some refiection, which will be extremely 
small unless the metric has short wavelength features.) The rest of the wavepacket 
will reach within a distance of order from the horizon, where it will undergo 
mode conversion and head back out, ending up in the constant velocity region as a 
combination of positive and negative free fall frequency wavepackets 1^+ ciiid The 
mode conversion process will be explained in the next section. 
Since the inner product is time independent we have 

{^out A) = {AJ) + 0'+ , 0) + ii'- , 0) , (44) 
or in terms of annihilation and creation operators, 

a(^out) = a{^r) + a{'ip+) - a^(^* ). (45) 

We assume that the state of the field at early times is the free-fall vacuum, |ff), which 
satisfies a{ip)\S) = for any ingoing positive free-fall frequency wavepacket ip- The 
particle creation in the packet ipout, characterized by the expectation value of the 



number operator, is thus given by the norm of tp^: 



(46) 



6. The Hawking effect and mode conversion 

In this section I will first describe the results of the particle production calculation 
in the Unruh model. These results will then be explained with the help of the WKB 
approximation. The resulting picture has two essential features: reversal of group 
velocity without reflection, and "mode conversion" from one branch of the dispersion 
relation to another. Interestingly, both these phenomena can occur for linear waves in 
inhomogeneous plasmasS^ii^and undoubtedly occur in many other settings as well. 
At the turn-around point the pure WKB approximation breaks down, and partial 
mode conversion from a positive free-fall frequency to a negative free-fall frequency 
wave takes place. This mode conversion gives rise to the Hawking effect. 

Unruh's computation^ starts with a purely outgoing, low wavevector, positive 
Killing frequency wavepacket in the asymptotic region where v{x) is constant (see 
top left graph in Fig. |^ for the wavepacket, top right graph for its power spectrum, 
taken from Ref. 0). Numerically integrating the partial differential equation (|3^) 
backwards in time, he finds that this wavepacket bounces off the horizon (middle left 
graph of Fig. |^) and moves back out to the asymptotic region with only very high 
wavevector components of both signs (bottom graphs of Fig. |^). The sign of the 
wavevector coincides with the sign of the free-fall frequency, so this 'final' wavepacket 
(which is actually the initial wavepacket going forward in time) possesses a negative 
norm component as seen from the power spectrum in the bottom right graph of Fig. 
0. Evaluating this norm Unruh obtains the occupation number (|^) for the outgoing 
packet. In all cases, the resulting occupation number agrees, to the accuracy of the 
computation, with the thermal prediction of Hawking. The computation contains alot 
more information than just this one number however, since a complete description of 
the initial wavepacket is obtained. Since the field equation is linear, one can simply 
fourier analyze the initial and final packets and read off the Bogoliubov coefficients 
for all wavevectors. Unruh did just this, and again found agreement with the ideal 
Hawking effect. 

These results have since been confirmed in two further calculations. In one 
approach!^, Brout, Massar, Parentani and Spindel (BMPS) introduced an analyti- 
cal approximation scheme in which the field equation near the horizon is solved by 
fourier transform. Using this technique, they showed that the spectrum of emitted 
radiation is thermal in leading order. In another approach, Corley and Jacobson0 
adopted the dispersion relation (|38[) and exploited stationarity of the metric to reduce 
the PDF to a fourth order ODE for each Killing frequency, which was then solved 
either numerically or, in some cases, analytically. This technique made it possible to 
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Fig. 7. A wavepacket bouncing off tfie horizon in tlie Unruh model, taken from Unruh's numerical 
integration results. The dispersion relation is (^7|), with ko = 512 and n = 1, and the Hawking 
temperature is T = 1.37. On the left is the wavepacket at three different times, with time going 
forward from bottom to top. The horizon is located at about x = 0.704, and the black hole lies to 
the right. On the right are power spectra in k of the initial and final wavepackets. 



go beyond leading order and compute the deviations from a perfectly thermal spec- 
trum. These deviations are relatively small at low frequencies where the Hawking 
radiation is copius, but are so large at larger frequencies that they eventually convert 
the exponential tail of the Planck spectrum to an increasing function. Nevertheless, 
the total flux associated with these higher frequencies is very small. 

6.1. Group velocity reversal and mode conversion 

The behavior of a wavepacket propagated back in time can be understood qual- 
itatively as follows. Let = e~*('^*~'^^\ substitute into the equation of motion, and 
drop terms containing derivatives in the free fall velocity. This is equivalent to using 
the WKB approximation, and yields the approximate position dependent dispersion 
relation 



This is just the dispersion relation in the local free-fall frame, since the free-fall 
frequency u' is related to the Killing frequency u by 



The position-dependent dispersion relation is useful for understanding the motion of 



who also found a Hamiltonian formulation for the wavepacket propagation 
using Hamilton- Jacobi theory. 

Graphs of both sides of equation (|47| ) are shown in figure ^ for F{k) given by 



{u-v{x)kY = F\k). 



(47) 



uj' = u — v{x)k. 



(48) 
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Fig. 8. Graphical solution of the position-dependent dispersion relation ( |47| ) with F{k) given by 
fTl). 



(|37| ) and for two different values of v. As x varies, the slope —v{x) (= of the 

straight line representing the left hand side of changes, but for a given solution 
the intercept uo is fixed since the Killing frequency is conserved. For a given x, 
the intersection points on the graph correspond to the possible wavevectors in this 
approximation. The coordinate velocity dx/dt of a wavepacket is the group velocity 
Vg = du/dk, which may also be expressed, using (|^), as Vg = du' / dk+v{x). Note that 
v'g = duj'/ dk is the group velocity of the packet in the free-fall frame, and corresponds 
to the slope of the curved line in figure 

Now assume the free fall velocity function asymptotes (at large positive x) to a 
value Vo satisfying —l<Vo<0, and consider a narrow wavepacket located far from 
the hole, centered about frequency u, and containing only k values around the smaller 
positive root (intersection point) shown in the figure. This is an outgoing wavepacket, 
since duo' /dk > \v\. Therefore, going backwards in time, the packet moves towards 
the hole. As x decreases \v{x)\ increases, so the slope of the straight line increases, 
until eventually the straight line becomes tangent to the dispersion curve. At this 
point Vg drops to zero. If uj is very small compared to /cq, then this stopping point Xt 
occurs when v{x) is very close to —1, that is, just barely outside the horizon. 

What happens at the stopping point? It was incorrectly suggested in Ref. that 
the wavepackets just asymptotically approach limiting position xt and wavevector kt. 
However, at the stopping point, the simple picture of a wavepacket peaked around 
a single wavevector has broken down. The wave packet has a significant spread in 
wavevectors which fail to satisfy the local dispersion relation. Furthermore, in this 
region there are nearby solutions to the dispersion relation. Just outside Xt, the 
straight line cuts the dispersion curve in two nearby points which approach each 
other as x approaches Xt- Simple WKB analysis is unable to say what happens. But 
a supplementary consideration gives the answer. As pointed out by Unruhlla, this is 
an unstable situation: for k slightly above kt the group velocity drops below zero (i.e. 
the comoving group velocity drops below the magnitude of the free-fall velocity) so, 
backwards in time, the wavepacket tends to move back away from the horizon. Once 
this begins to happen, k continues to increase as the wavepacket moves further away. 
Exactly this behavior was found in Unruh's numerical solutionS to the PDE. 

This is not the end of the story however. There is yet another "nearby" solution 
to the dispersion relation, on the negative wavevector branch, which is excited as 
well. This process is an example of "mode conversion" 110. The mode conversion is 
strong if the wavelength and phase velocity of one mode is close to that of another. 
It appears that this condition for strong conversion is not met, since the negative k 
solution kt- to (^7]) at xt is not very close in magnitude to kt (see Fig. ||.) However, 
one must remember that the wavepacket here has a significant spread in wavevectors 
which fail to satisfy the local dispersion relation. Thus the wavepacket can pick up 
other modes that fulfill the conversion criterion. That the wavepacket must pick 
up a negative wavevector piece is clear in the case of the ordinary wave equation. 



since there causality implies that the wavepacket strictly vanishes inside the horizon, 
which cannot be accomplished with only positive wavevectors. In the Unruh model, 
the wavepacket satisfies approximately the ordinary wave equation until it gets very 
close to the horizon, so one should expect that a similar amount of mode conversion 
takes place. Evidently the negative wavevector mode does couple in strongly for 
sufficiently small cj, as shown both by Unruh's solution of the PDE and by the ODE 
methods applied by BMPS0 and ourselveJ^ll. The "converted" , negative wavevector, 
wavepacket also has a negative group velocity and so also moves, backwards in time, 
away from the hole. The end result thus consists of two wavepackets, one constructed 
of large positive k wavevectors and the other of large negative k wavevectors, both 
propagating away from the hole and reaching the asymptotically fiat (constant free 
fall velocity) region. Note that the negative wavevector component will consist of 
wavevectors of magnitude slightly larger than the positive wavevector component. 
This is clear from the position of the intersection points in Fig. ^, and is borne out in 
Unruh's computation (Fig. ^ as well. The number of created particles in the final, 
late time, wavepacket is given by the norm of the negative wavevector part of the 
initial, early time wavepacket. 

7. The stationarity puzzle 

Particle production via the Hawking effect in the case of a black hole that forms 
in a collapse process is ordinarily understood to depend critically on the fact that 
the metric is not static for all times. Although the radiating modes come in from 
infinity and go back out to infinity, their Killing frequency is not conserved because 
they propagate through the time dependent part of the spacetime. The way the 
usual Hawking effect transpires in a strictly stationary spacetime is that the outgoing 
wavepackets are traced backwards to parts that do not make it back out to infinity, 
but rather cross the white hole horizon, at which point the Unruh boundary condition 
on the quantum state is imposed. The piece of the wavepackets that scatters off the 
curvature and does make it back out to infinity is not associated with particle creation. 
If, as in the Unruh model, the entire wavepacket turns around and goes back out, 
remaining for all time in the stationary region of the spacetime, then it would seem 
that there can be no particle creation at all. So how does the Unruh model yield a 
nontrivial Hawking fiux? 

The first answer is that the wavepackets have not been followed all the way out 
to infinity. Another way to put it is that we have taken the free fall frame, in which 
the boundary condition is imposed, to be moving towards the black hole at infinity, 
rather than coinciding with the rest frame of the black hole at infinity. (Technically, 
this corresponds to the fact that the metric function v{x) does not vanish at infinity.) 
But what happens if we take v{x) to vanish at infinity, and continue to follow the 
wavepackets backwards in time, not stopping to impose the quantum state boundary 



condition until the packets reach infinity — or do something else? Do the wavepackets 
propagating backwards in time ever reach infinity? Is there any Hawking radiation? 

What happens in the Unruh model with the dispersion relation ( p7| ) is that the 
magnitude of the wavevector grows without bound as the wavepacket moves outward 
where v{x) is falling to zero. Thus, even though the difference between the free- 
fall and Killing frames is going to zero, the wavevector is diverging in such a way 
that the wavepacket always maintains a negative free-fall frequency part of the same, 
negative, norm. Thus the Hawking effect indeed occurs. From this analysis we see 
that the Unruh model, while it entails a strict cutoff in free-fall frequency, involves 
in an essential way arbitrarily high wavevectors, i.e., arbitrarily short wavelengths. 
Insofar as we wish to explore the consequences of a fundamental short distance cutoff 
on the Hawking effect, this is an unsatisfactory feature of the model. The outgoing 
modes emerging from the black hole region still arise from arbitrarily short wavelength 
modes. Using the dispersion relation (^) on the other hand, the wavevectors are 
bounded by /cq. In this case something very strange seems to happen0, which remains 
to be understood: the positive norm piece of the wavepacket goes backwards in 
time out to infinity at superluminal group velocity, and the velocity of the negative 
norm piece goes past positive infinity to negative infinity, which seems to result in a 
propagation "back to the future" . 

In search of a more realistic model, it is interesting to go back and consider the 
behavior of a wavepacket propagated backwards in time using the dispersion relation 
of liquid helium-4 (see Fig. ||). Up to the reversal of group velocity outside the 
horizon, the behavior is the same as for the Unruh model. After that (i.e. before 
that) the packet will go over the first maximum of the dispersion curve, at which 
point its comoving group velocity changes sign, which only pushes it away from the 
horizon even faster. Eventually, however, it approaches another turn around point, 
near the roton minimum, where the free-fall frequency line becomes tangent once 
again to the dispersion curve. It seems reasonable to suppose that what happens 
here is another reversal of direction and further mode conversion: the wavepacket 
continues along the dispersion curve, falling back towards the horizon, until it runs 
off the end of the quasi-particle spectrum, where it is presumably unstable and has 
a wavelength on the order of the interatomic spacing. Or, perhaps, the instability 
sets in even before this point. It thus appears that there is no way to analyze the 
helium model fully at the level of quasi-particle field theory. Rather, the many particle 
dynamics must be directly confronted. It would be very interesting, though perhaps 
very difficult, to analyze this many body problem. 
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